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ABSTRACT 

We study a class of non-local, purely gravitational models which have the 
correct structure to reproduce the leading infrared logarithms of quantum 
gravitational back-reaction during the inflationary regime. These models end 
inflation in a distinctive phase of oscillations with slight and short violations 
of the weak energy condition and should, when coupled to matter, lead to 
rapid reheating. By elaborating this class of models we exhibit one that has 
the same behaviour during inflation, goes quiescent until the onset of matter 
domination, and induces a small, positive cosmological constant of about the 
right size thereafter. 
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• Introduction: During tlie inflationary era infrared gravitons are pro- 
duced out of the vacuum because of the accelerated expansion of spacetime. 
The interaction stress among the gravitons produced - an inherently non- 
local effect - can lead to a non-trivial quantum gravitational back-reaction 
on inflation pTj. In a previous paper |2] we proposed a phenomenological 
model which can provide evolution beyond perturbation theory. In one sen- 
tence, we constructed an effective conserved stress-energy tensor T^u[g\ which 
modifles the gravitational equations of motion: 

G^u = R^iu - \^9iivR = -^gfiu + 8nGT^i,[g] . (1) 

and which, we hope, contains the most cosmologically significant part of the 
full effective quantum gravitational equations. 

Our physical ansatz consisted of parametrizing T^j,^ [g] as a "perfect fluid" : 

Ti,4g] = {p + p)Uf,u^ + pQf,^ , (2) 

with the gravitationally induced pressure given as the following functional of 
the metric tensor: 

p[g]{x) = A'f[-eX]{x) , ^ ^ , (3) 

where the function / grows without bound and satisfles: 

f[-eX] = -eX + 0{e') , (4) 
and where the scalar d'Alembertian: 

□ = -^d,{g'^^^d,) , (5) 

is deflned with retarded boundary conditions. The induced energy density 
p[g] and 4- velocity u^[g] were determined, up to their initial value data, from 
stress-energy conservation: 

D^T^. = . (6) 

^Hellenic indices take on spacetime values while Latin indices take on space values. 
Our metric tensor g^^^, has spacelike signature and our curvature tensor equals: R°'p^^ = 
r'l/s,;. + r"^p - (m ^ v)- The initial Hubble constant is 3iJ^ = A. We restrict 
our analysis to scales M = (A/87rG')i below the Planck mass Mpi = G^i so that the 
dimensionless coupling constant e = GA of the theory is small. 



1 



The 4-velocity was chosen to be timehke and normahzed: 



g>'''u^u, = -l =^ u>'u^-u = ^ ■ (7) 

We cannot exhibit u^[g\{x) and p[g\{x) as exphcit functionals of the metric 
the way expression ([3]) does for the pressure, but it is straightforward to solve 
the conservation equations for a homogeneous and isotropic background El and 
even for perturbations [3j. For the analysis that will follow we require only 
the background results: 



u^it) = -a{t)6' , p{t) = -Pit) + / dt'Pit' 







a{t) 

Hence the equation of state of our ansatz for the induced stress-energy is: 



(8) 



pit) 

-^'^ - W) 



_ :^ pit') a\t') 
Pit) a\t) 



(9) 



The homogeneous and isotropic evolution of this model - using a com- 
bination of numerical and analytical methods - revealed the following basic 
features El : 

01. After the onset and during the era of inflation, the source X{t) grows 
while the curvature scalar R{t) and Hubble parameter H{t) decrease. 

02. Inflationary evolution dominates roughly until we reach a critical point 
Xcr defined by: 

1 -SyrGA/i-GAXe,] = . (10) 

03. The epoch of inflation ends close to but before the universe evolves to 
the critical time. This is most directly seen from the deceleration parameter 
since initially q{t = 0) = — 1 while at criticality q{t = tcr) = +|. 

04. Oscillations in R{t) become significant as we approach the end of infla- 
tion; they are centered around R = 0, their frequency equals: 



u = GAHoJ72nn, , (11) 



^The line element in co-moving coordinates is ds^ = —dP + a'^{t) dx-dx. In terms of the 
scale factor a, the Hubble parameter equals H{t) = da^^ and the deceleration parameter 
equals q{t) ~ —acid^^. 

^In [2] , our analytical results were obtained for any function / satisfying ([4]) and growing 
without bound, our numerical results for the choice: f{x) = exp(a;) — 1. 
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and their envelope is linearly falling with time. 

05. During the oscillations era, although there is net expansion, the oscil- 
lations of H{t) take it to small negative values for small time intervals - a 
feature conducive to rapid reheating; those of H{t) take it to positive val- 
ues for about half the time; and, those of a{t) are centered around a linear 
increase with time. 

Since we shall be concerned with the various post-inflationary phases to 
which the universe evolved, and since all these phases are characterized by 
constant e = —HH~^ = 1 + g , it is convenient to provide expressions for 
some basic quantities describing such spacetimes. In particular, in terms of 
their initial values, the time evolution of the scale factor, Hubble parameter 
and Ricci scalar are, respectively: 

^^^^ ^ a{t) ^ aiJl + eHir,{t-Un)Y , (12) 



^ ^ 1 + eHin{t - tin) ' ^ ^ 

R{t) = 6H + 12H^ = 6{2-£)H\t) . (14) 

We shall also need the temporal component of the Ricci tensor: 

i?oo = +SqH^ = -3(1 -e)H^ . (15) 

Moreover, the first two time derivatives of the curvature scalar are: 

R = -l2e{2-e)H^ , (16) 

R = +36£2 (2 -e)H^ , (17) 

so that: 

UR = -R-3HR = +36e {1 - e){2 - e) . (18) 

Notice that R vanishes for radiation (e = 2) and HR for radiation as well as 
inflation (e = 0) . 



• Two Problems: The homogeneous and isotropic evolution described 
by the simple model briefly reviewed in the introduction does not give a 
completely satisfactory end to inflation. The oscillations that occur after the 
end of inflation are not a problem, but the average expansion a{t) oc t is 
unacceptably rapid. At that rate there would be no reheating and the late 
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time universe would be cold and empty. Nonetheless, the same is true for 
scalar-driven inflation if one ignores the possibility for energy to flow from 
the inflaton into ordinary matter. We believe that energy will flow from the 
gravitational sector of our model into ordinary matter to create a radiation 
dominated universe, just as it is thought to do for scalar-driven inflation. 

An amazing possibility arises if this process can be shown to occur: our 
quantum gravitational correction cancels the bare cosmological constant and 
then becomes dormant during the epoch of radiation domination. To see this, 
suppose the deceleration parameter has the pure radiation value of q{t) = +1 
for times t > tr- This case corresponds to e = 2 and (11211141) give: 



+1 



a{t) 

m 
m 



Or 1 + 2Ho{t - tr 

l + 2Hr{t-tr) ' 
, 



(19) 

(20) 
(21) 



Our simple source X{t) obeys the differential equation CX = i?, so for t > tr 
it must be a linear combination of its two homogeneous solutions: 



Vt > tr 



□X = 



X{t) = Xr + Xr dt' 





Or 


1^ dt' 




_a{t')_ 



^r — 



Hr Jl + 2Hr{t - tr) 



The only solution consistent with g = +1 is: 



Xr 



Xr 



Xr 







(22) 



(23) 



Note from equation ([3]) that X{t) becoming constant implies that the induced 
pressure pit) is also constant. One can see from equations ([8ll9]) that the 
energy density also becomes constant provided the approach to radiation 
domination results in the following condition: 







dt p{t) a^{t) 







(24) 



^For all power laws, such as radiation, the pressure and density fall like t^^ . Therefore, 
the second of the two homogeneous solutions cannot be present, implying that Xr — 0, 
since its time dependence is t^^ and could not sustain radiation. Neither can the first 
homogeneous solution sustain radiation unless the constant Xr eliminates the cosmological 
constant iH^ and, hence, equals Xcr- 
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In fact this is also the condition for achieving the fixed point solution 
and it results in our induced stress-energy becoming perfect vacuum energy 
which completely cancels the bare cosmological constant A. 

Having X{t) approach Xcr within the context of a hot, radiation dom- 
inated universe would be a great success for our model, but the eventual 
transition to matter domination poses enormous problems. The onset of 
matter domination is really a gradual process but let us simplify the expo- 
sition by considering a sudden change from q = +1 to q = +^ at some time 
tm ^ U- During this matter dominated epoch, for which e = |, expressions 
(fiaiTl become: 



+ - 



1 



a{t) 

m 
m 



1 + -j^Hm{t — U, 



1 -I- ^Hmit — U. 



3H^ 



1 -|- 2 Hffi (t 



(25) 
(26) 

(27) 



where Hm and H{tm) and a{tm), respectively, computed from the 

radiation dominated geometry fll9H20p . The resulting change in the source 
X{t) is: 



q = + 



2 

AX{t) 



X{t) - X, 



In 



1 + -Hmit — tr, 



0(1) . (28) 



The fact that matter domination causes X{t) to move off its critical value X^r 
means that our induced stress-energy no longer behaves like perfect vacuum 
energy. That would pose no problem for cosmology provided that: 

(i) the residual effect from the perturbation AX(t) results in a pressure and 
energy density which are initially small compared to the energy density of 
matter, 

(ii) the residual effect eventually acts to induce a phase of acceleration after 
a suitably long time lag. 

To understand what is wrong with the change fl28|) caused by matter 
domination, it is useful to recall our ansatz for the quantum gravitationally 
induced pressure: 

p[g]{x) = A\f[-GAX]ix) . (29) 
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In the context of this ansatz there are two major problems with 

01. The sign problem. It derives from the function f{x) in (!29l) be- 
ing monotonically increasing and unbounded. Hence, pushing X{t) below 
Xcr <^ results in positive total pressure, whereas observation implies nega- 
tive pressure during the current epoch [H [5]. Note that we cannot alter this 
feature of f{x) without sacrificing the very desirable ability of the model to 
cancel an arbitrary bare cosmo logical constant. 

02. The magnitude problem. In one sentence, the magnitude of the total 
pressure produced by fl28|) is vastly too large. The problem arises from the 
factors of the bare cosmological constant A in our ansatz (121?]) . The total 
pressure ptot is the sum of the classical contribution and our ansatz fl29]) : 

Ptot = - ^ |l - 8^^^ f[-GA {X,r + AX) ] I (30) 

^ - ^ X (GA)2 AX . (31) 

Comparing with the currently observed value Pnow of the pressure: 

3 



SttG 



Hnow ) (32) 



gives: 



, /; AX ^ 10«6 X X AX , (33) 

where we have assumed Hq ~ 10^'^ GeV and i/now ~ 10^^'^ eV. The deriva- 
tive /^^ is unity for the linear model and of order (GA)^^ ~ 10^^ for the 
exponential model, so we expect /^^ to be at least of order one and possibly 
much greater. 

There is no way of addressing either problem without generalizing our 
ansatz fl29|) for the pressure. This necessarily takes us away from what can 
be motivated by explicit computation during the de Sitter regime. 

• Decreasing the Magnitude: The magnitude problem arises because 
the constant A in ( |29|) is about the square of the inflationary Hubble parame- 
ter rather than its late time descendant that could be 55 orders of magnitude 
smaller. Solving the problem entails replacing one of these factors of A by 



6 



some dynamical scalar that changes as time evolves in a way that also pre- 
serves the original relaxation mechanism. The latter requirement rules out 
any tampering with the factor of that multiplies /[— GAX] in f l29|) . It is 
safer to make the factor of A in the argument of the function / dynamical 
and move it to the right of the D"^: 

-GKX = "^(^""x) ^ -G^^"^[9] ' (34) 

where S is an appropriate scalar quantity. The idea is for ( l34|) to approach 
— GA X Xcr during inflation and then freeze in to this value during radiation 
domination, throughout which the scalar falls off so that subsequent evolution 
is driven by an acceptably small source. There are many possibilities for the 
scalar 5" in (IMj) . which we shall always normalize so as to make the ratio 
SA~^ give unity for de Sitter spacetime. 

One might think S = can work, but numerical simulations show that 
it responds too quickly to the slowing geometry. Instead of inflation ending, 
the scale factor approaches an accelerating, power law expansion - that is, 
a{t) oc with s > 1. It is instructive to understand why this happens 
by perturbing the pressure around the critical value which cancels the bare 
cosmological constant. The relevant FRW equation of motion becomes: 

-2H-3H^ = -A + 8nGp (35) 

= -A I 1 - 8nGA f [ - GA{Xcr + AY) ] | (36) 

^ -87T{GAy f'^^ AxAY. (37) 

When S = ^R, the function AY obeys: 

1 R"^ 
S = -R =^ HAY = — . (38) 
4 4A ^ ^ 

Then, for a{t) oc we have: 

a{t) oc f ^ R = - ^ ' 



Yi 3s\2s-lf , , 
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where Yi is an integration constant. For s > 1 the homogeneous solution 
proportional to Yi falls off faster than the inhomogeneous solution. Neglect- 
ing the homogeneous solution and substituting into fl371) gives an algebraic 
equation for the power s, which is easy to solve for GA <^ 1: 



TTT^ = 167r(G'A) 



s[s- 2> 



1 



167r(G'A)2/^, 



> 1 , GA<1 . (40) 



A more non-local scalar, which responds less rapidly to the slowing ge- 
ometry, can be formed from derivatives of X: 

S _ -g^'^ d,X[g] d^X[g] X' _ 1 

A" fA ^ fA ' M9] = ^R- (41) 

This gives an end to inflation for relatively large values of GA but still goes 
over to a power law with the asymptotic form fHOj) for small GA. The reason 
is obvious from the form of X for FRW: 

X — ^ f dt' o!^it')Rit') . (42) 
a'^(t) Jo 

If R simply vanished after some time, then X would decay subsequently like 
a~^{t). Hence the lag to changes in R is of the order of a few Hubble times 
and inflation will only be ended for unrealistically large values of GA such 
that the critical time tcr falls within this period. 

We also tried scalars formed by using inverses of other operators like the 
conformal d'Alembertian: 

1 „ 1 d d , , , 

n,^^n--R ^ («) 

and the Paneitz operator [6l [7] : 

Dp ^ □^ + 2Z^Ji?--Vi?)D. -> -A.-i-i4.- (44) 



3 ^ a? dt dt dt dt 

However, the result was always power law expansion like f l40|) for GA ^ 1. 

To avoid the magnitude problem and still end inflation requires that we 
evaluate the dynamical scalar far back in the past. One way of achieving this 
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is to use the integral curves XtAd]^^) ^ timelike 4- velocity field We 
can construct such a 4- velocity field by taking the gradient of the invariant 
4- volume V[(7](x) of the past light-cone from the point = {t,x) back to 
our initial value surface at x''^ = {t,x') [8]: 

V = f^dt' j d^x' ^-g{x') 6(^- a{x]x')) . (45) 

Here a{x] x') is | times the square of the geodesic length from x^ to x'^ 
[9]. Because the volume of the past light-cone grows as one evolves into the 
future, its gradient is guaranteed to be timelike in any geometry. We define 
V^^[(7](a;) as the normalized gradient of V[5'](x): 

V^[g]{x) = , ^ . (46) 



Now construct the integral curves x'^[(7](t, x) - as a functional of the metric 
and an ordinary function of the parameter r and a coordinate point x^ - so 
that they obey the conditions: 

^ = -V^{x) and x^(0,a;) = . (47) 

Our physical requirement on the scalar S is to evaluate it at any time t* such 
that infiation is still dominant. Suppose T*[g\{x) corresponds to this time as 
we follow the integral curves from x^ back to the initial value surface. Then, 
we might define the scalar S at point x^ in expression ( l3l|) to equal: 

S{x) = ^R[x^{t*;x)] . (48) 

The modification we have outlined would not change the initial value 
problem, is defined for any geometry, and seems to be invariant. It also 
solves the magnitude problem while preserving the general features of the 
way the original ansatz ( 129|) ends infiation. For an FRW geometry and for 
the value of T*[g]{x) associated with, e.g. ^ of the time from x^ back to the 
initial value surface, the functional y[(7](x) becomes: 

i-blW ^ -J^,fJ^J^',,"aV')Rin.^. (49) 
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Figure 1 shows the Hubble parameter versus time for this model with GA = 
Any of the less non-local models we tried evolved to the power law 
solution f HOj) for this small a value of GA. The duration of the plot corre- 
sponds to 350 initial Hubble times. One can see that the period of oscillation 
lengthens towards the end because j^Rijot) decreases, but H{t) still oscil- 
lates with decreasing amplitude and the oscillations still drop below zero. For 
the enormously smaller values of GA relevant to primordial inflation -^Rijot) 
would be effectively constant during the oscillatory phase, so the period of 
oscillation would also be constant, just as in the simple model (!29|) . 

To recapitulate, we solve the magnitude problem by changing the con- 
stant A in GA D~^R into a factor of the dynamical quantity R, standing to 
the right of the inverse d'Alembertian. When the eventual onset of matter 
domination perturbs our induced stress-energy off its radiation dominated 
fixed point, the change will be small - rather than overwhelming - because 
the multiplicative factor of A has been replaced by the vastly smaller value 
of R at nearly the time of matter domination. However, evaluating this ad- 
ditional factor of R at the same time as the original factor of R makes the 
source "turn off" too quickly to end primordial inflation. Our solution is 
to instead evaluate the additional factor of at a time T*[5f](x) which is 
early enough that the Ricci scalar of primordial inflation is still nearly 4A. 
Although we set r* to be one tenth of the current time for the sake of defi- 
niteness, the only two requirements on it are: 

(i) Near the end of primordial inflation the time r* must be early enough so 
that we still have R{t*) ^ 4A. 

(ii) After the onset of matter domination r* must be late enough so that 
R{t*) is comparable to the square of the Hubble parameter at that time. 
There are an enormous number of functionals r[(7](a:) with these two proper- 
ties and any of them would suffice. 



• Changing the Sign: We turn now to the sign problem. It arises 
because the Ricci scalar ( !T^ is positive during both inflation (e = 0) and 
matter domination (e = |): 

q = -I =^ R = +12H^ , (50) 

q = +1 =^ R = +SH^ . (51) 

What we need is a source which changes sign from inflation to matter domi- 
nation, and is still zero (or very small) during radiation domination. There 
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are again many possibilities but a simple one that works is to change the 
source from the form motivated by 



Y[g]{x) = - 



R{x) X 



R{x{r*,x)) 
4A 



to the following form: 



Ya[9]{x) = -^^[R{x)R{x{r\x)) + anR{x) 



(52) 



(53) 



The term proportional to a vanishes for both de Sitter inflation (e = 0) and 
pure radiation domination {e = 2), so it should make little difference until 
the onset of matter domination. For t ^ tm and for the choice of T*[5r](a;) 
already discussed, the various dynamical scalars in (!53|) give (e = |): □ 

t > =^ R(t) = +3H^{t) , R{^t) = +100 X 3H^{t) , (54) 

27 

nR{t) = -^H\t) . (55) 

where we also used f|T8l) . There exists a reasonable range of a for which the 
sign gets reversed. 



• An Even Simpler Ansatz: Although expression ( l53l) resolves the 
sign problem, it introduces a new parameter a to the theory and, more 
importantly, it adds initial value data since the second term contains 
fourth order derivatives of the metric. |^ 

Both issues can be avoided by adopting the following geometrically mo- 
tivated form for the source: 



V^ix{T\x)) V^(x{r\x)) RJx{r\x) 



(56) 



which - for an FRW geometry and for the value of r* [g\ [x) associated with, 
e.g. ^ of the time from x^ back to the initial value surface - becomes: 



Ys[9]{x) 



dt' 



a3(t') 



dt" a%t") R{t") X 



R 



A 



(57) 



^Because during matter domination H{t) = and R{t) = 3H^{t), it trivially follows 
that: Ri^t) = 10^R{t). 

^Unlike the first term jR^ which has the square of second order derivatives acting on 
the metric. 
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where we have used that for these particular spacetimes: 



FRW 



V^[g]{x) 



(58) 



• Late Time Acceleration: The ansatz (156|) and its FRW form (|5j 
has essentially the same time evolution with that of the simple model ( 129|) 
during inflation, vanishes during radiation domination, and the transition 
to matter domination can trigger a slow, very small decrease in p[g\. This, 
in turn, can uncover a tiny portion of the bare cosmological constant and 
eventually cause a phase of acceleration. 

To see this behaviour quantitatively, note that: 



R{t) X Rooi^t) = X 300i/^ = 
= Si/^ X 300//^ = 
= 3i^2 X 150/f2 = 



900i/^ 
450i^^ 



where we have used (11411151) . 

In view of (16011611) . to compute the source Y 



10 tr <t <tm 
10 tfjl *\ t <C tjiow 



(59) 
(60) 
(61) 



x) given by ([56 



Ys[gm 



1 1 
An 



(62) 



we must evaluate the action of the inverse differential operator on H'^{t)] 
we might as well do so for any constant e and for a generic lower limit of 
integration which we shall denote by T: 



□ 



3(1-5 
1 

'3(1-8 



" 27 ^ + 27 ^ + 3^ a^{t)H{t) 

H^{T)\ . 



(63) 



3-e 

This is the full answer for any e. Of the four terms present in ( l63l) . the 

2' 



first and third are sub-dominant at late times for matter domination (e = |) 
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because of their t ^ and t ^ behaviour, respectively. The remaining two 
terms in fl63p are constant and give: 



t > U 
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(64) 



Taking into account fl^T]) and setting the lower limit of integration equal to 
T = tm, the relevant part of the source (!62l) becomes: 



t > 



200 H, 



(65) 



where Hm = H{tm)- 

In terms of the new source l5[g'](t), the total pressure fl30|l becomes: 



A 



SttG 
A 

8^ 



[l-^TvGkf[-GkYs\} 
\ 1 - SttGA / [ - GA {X^r + AYs) 



t > =^ p,^, = _ __ ^ 1 _ SnGA f[-GA Yg] \ (66) 

(67) 

(68) 
(69) 

where in the last step we used fl65|) . The ratio of the predicted total pressure 
( 169|) to the current total pressure ( l32l) is: 



200 GA' f„ Hi , 



t > U 



Ptot 
Pnow 



^ 8vr(GA)^ X /; X 10^° 



For the exponential model analyzed in [2j, we have: 

f(x) = ^ 
SO that the pressure ratio (fTTI) takes the form: 



1 



SttGA 



~ 3 ^ ^ 



(70) 

(71) 

(72) 
(73) 



Physical values of the dimensionless coupling constant GA = M'^Mpf easily 
balance the factor 10^^ so that we can achieve equality of the predicted total 
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pressure ptot to the observed total pressure and, hence, account for the current 
acceleration of the universe. 

To recapitulate, the onset of matter domination causes our improved 
ansatz for the quantum stress-energy to move slightly off the point at which 
it exactly cancels the bare cosmo logical constant. We have shown that the 
shift approaches a small, positive cosmological constant which could easily 
have the right magnitude to induce the current phase of cosmological acceler- 
ation. For most of the matter dominated phase this solution is self-consistent 
because the residual cosmological constant is too small to alter the expan- 
sion history until very late. Furthermore, the eventual onset of acceleration 
induces a new change in y5[5(](t), this time in the sense of diminishing the 
residual cosmological constant. However, the time scale for that to occur is 
the current Hubble scale, so one would not expect much effect as of yet. 

• Epilogue: We have presented an improved ansatz for the most cos- 
mologically significant part of the effective field equations of quantum gravity 
with a positive cosmological constant. This improvement of the simple ansatz 
analyzed in [2] and reviewed in the introduction herein, allows us to expect 
a reasonable post-inflationary time evolution into the present. 

Prominent among the list of topics for future work is perturbations [3]. 
We need to show that the dynamical scalar mode of our model dumps the 
energy of oscillations into matter to reheat the universe in a natural way. 
If this happens, the improved model presented here can make the quantum 
gravity sector go quiescent during a long epoch of conventional radiation 
domination. The subsequent transition to matter domination might even 
give rise to something like the current phase of acceleration without severe 
fine tuning. 

Furthermore, we must derive and solve the equation for scalar perturba- 
tions, at least enough to compute the scalar power spectrum. In contradis- 
tinction, the equation for tensor perturbations remains unchanged and we 
need only use the expansion history a{t) predicted by our model in order to 
compute the tensor power spectrum. Finally, one also needs that there be 
no long-range scalar force at late times. 
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Figure 1: The Hubble parameter (in units of Hq) versus time (in units of j^qqHq) 
for p = A^f[-GAY] with GA = ^ and Y given in 
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